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ABSTRACT 

The angular power spectra of the infrared maps obtained by the DIRBE (Diffuse 
InfraRed Background Experiment) instrument on the COBEQ satellite have been 
obtained by two methods: the Hauser-Peebles method previously applied to the DMR 
maps, and by Fourier transforming portions of the all-sky maps projected onto a plane. 
The two methods give consistent results, and the power spectrum of the high-latitude 
dust emission is Ci oc £~ s in the range 2 < £ < 300. 

Subject headings: ISM: structure; dust; cosmic microwave background 

1. Introduction 

The angular power spectrum of the sky brightness at millimeter wavelengths can be used to 
determine the parameters of cosmological models, such as H Q , f2 Q , £Ib and A (Jungman et al. 
1996), but our ability to measure this power spectrum will depend on the angular power spectrum 
of foreground emission from the Milky Way. Kogut et al. (1996) have found that free- free emission, 
which will be the dominant foreground contamination, is better correlated with far infrared dust 
emission than with synchrotron emission. The angular power spectrum of the dust emission also 
determines the beam size dependence of the confusion limit to the sensitivity of point source 
surveys. Gautier et al. (1992) found a power spectrum P(k) oc k~ 3 for the IRAS 100 /j,m maps, 
in a range of k corresponding to 200 < £ < 3000, and Kogut et al. (1996) found Cg oc £~ 3 for the 
DIRBE (Boggess et al. 1992) 240 fim map in the 2 < £ < 30. But the DIRBE instrument has an 
angular resolution of 0.7°, which allows one to measure the power spectrum to angular frequencies 
much higher than £ = 30, and DIRBE has 10 bands (240, 140, 100, 60, 25, 12, 5, 3.5, 2.2 and 
1.6 /im) which allow one to measure the frequency dependence of the angular power spectrum. 
In this note I use a simple method to allow for the effect of detector noise, which dominates the 
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100 < 1 portion of the DIRBE 240 /um power spectrum, and show that the 60-240 /um DIRBE 
maps follow a Cg oc i?~ 3 spectrum over the range 2 < £ < 300. Since this spectrum is more than 
a full power of i steeper than the predicted spectrum of the cosmic microwave background in 
inflationary dark-matter dominated models, observations of the "Doppler" peak at I ~ 200 should 
be much less affected by Galaxy than the COBE DMR maps. 

2. Method 

The power spectrum of a process is a quadratic function of the data, so noise will bias the 
power spectrum even if it is not correlated with the signals. The technique used by Wright et 
al. (1994) to correct for this effect in the COBE DMR maps was to use the two sides of the 
DMR instrument, the A and B sides, to make sum and difference maps, S = (A + B)/2 and 
D = (A — B)/2. The noise power spectrum canceled out when the difference in the power spectra of 
the S and D maps was taken to yield the the power spectrum of the sky. The DIRBE experiment 
does not have a comparable duplication of detectors, but I have computed maps with independent 
detector noises by forming sums of the even weeks and the odd weeks of the data. The DIRBE 
data have been averaged into separate files for each of the 41 weeks of data taking to allow for the 
removal of the time varying emission from interplanetary dust. This time variation is due to the 
motion of the Earth around the Sun, which means that a given spot on the sky is viewed through 
different parts of the interplanetary dust cloud, and the time variation is used to determine the 
properties of the interplanetary dust. Once the zodiacal light is removed (see Appendix |A|), the 
weeks are summed into all-sky, all-data maps. By summing the even and odd weeks, I make E 
and O maps, and construct the sum and difference maps as S = (E + 0)/2 and D = {E — 0)/2. 
Figure [l] shows the sum and difference maps for 100 fxm DIRBE maps, while Figure || shows the 
240 fim maps. The pass 2B DIRBE weekly maps, available since 1994 on the NSSDC, have been 
used in this paper. Figures [l] and || have been histogram equalized, so the values on the scale 
bar correspond to the th , 20 th , A0 th , 60 th , 80 th and 100 th percentiles of the map. The difference 
map at 240 //m is clearly noise dominated, with the middle 60% of the histogram spanning 1.85 
MJy/sr. A Gaussian distribution with a standard deviation of 1.1 MJy/sr would have the same 
middle 60% span, but the actual distribution of values is not Gaussian due to variable coverage 
and to errors caused by large sky gradients near the galactic plane. This detector noise is a very 
significant contribution to the power spectrum of the sum map. The difference map at 100 /iin 
is dominated by a systematic error in the zodiacal light removal caused by an inadequate model 
for the interplanetary dust cloud, but the overall middle 60% span of 100 kJy/sr is very small 
compared to the sum map. The detector noise in the 100 /im channel on DIRBE is 100 times 
smaller than the detector noise in the 240 /mi channel (Boggess et al. 1992) but the middle 60% 
span of the difference map is only 19 times smaller, showing that the 100 /xm difference map is 
not detector noise dominated. Since the zodiacal light model error has a period of two weeks due 
to the alternating signs given to successive weeks in the difference map, the dominant spherical 
harmonic in this difference map is I = 26. The sum map will not have this I = 26 power because 
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BAND 8 Sum 




0.2880 1.9711 3.2075 5.4438 12.032 8556.6 



BAND 8 Difference 




-213.9 -.0430 -.0075 0.0187 0.0570 178.67 



Fig. 1. — Top: 100 /xm map from the sum of the DIRBE even and odd week sum maps. The North 
Galactic Pole is the center of the left circle, while the South Galactic Pole is the center of the right 
circle. Bottom: Difference of the even and odd sum maps. 



- 4 - 

BAND 10 Sum 




-3.481 2.6945 4.7228 8.2842 20.793 7056.0 



BAND 10 Difference 




-122.4 -.9443 -.2891 0.2586 0.9087 181.42 



Fig. 2. — Top: 240 fim map from the sum of the DIRBE even and odd week sum maps. Bottom: 
Difference of the even and odd sum maps. 
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Fig. 3. — Power spectra of the DIRBE 240 //m map after 4 different galactic plane cuts: the 
"Designer Cut", and \b\ > 30°, 40° and 49°. The lines are best fits with an assumed slope. 

it does not apply alternating signs to successive weeks, so this method of correcting for noise bias 
will not work for maps where the zodiacal light errors are comparable to or larger than the true 
signal. This limits my analysis to the A > 60 fim DIRBE maps. The 25 and 12 DIRBE maps 
are so dominated by the interplanetary dust emission that the residuals due to inadequate models, 
while only a few percent of the peak emission, are still larger than the high galactic latitude 
emission from interstellar dust. 

3. Variation with latitude cutoff 

The power spectrum of the DIRBE maps is a strong function of the galactic latitude cutoff. 
This indicates that the galactic emission is not a stationary Gaussian random process, even 
after the very bright emission from the Galactic plane has been cut out. Figure || shows the 
Hauser- Peebles (1970) estimate of the power spectrum of the DIRBE 240 //m map with various 
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cutoff latitudes. The lowest latitude cutoff is labeled DC for "Designer Cut" - this is the same 
part of the sky used for the analysis of the 4 year DMR maps (Bennett et al. 1996). This is 
basically a \b\ > 20° cut with some further exclusions in Taurus-Orion and Scorpius-Ophiucus, and 
uses 63% of the sky which is the area given by a straight \b\ > 22° cut. The power spectrum in 
this "Designer Cut" region is very similar in magnitude to the \b\ > 30° cut, even though the 63% 
sky coverage of the "Designer Cut" is quite a bit larger than the 50% sky coverage in the |fe| > 30° 
cut. However, the \b\ > 40° and especially the \b\ > 48.6° cuts give power spectra that are lower 
still but at a big cost in sky coverage (36% and 25% coverage respectively). 



4. High frequency power spectrum 

The Wright et al. (1994) method of finding power spectra is very slow for large £'s, so an 
alternative approach using the fast Fourier transform has been used for £ > 9. The equal area 
"polar eyeball" projection used in Figures [FJ and § is used to map one galactic hemisphere onto a 
square with 512 x 512 pixels. This map is windowed with the function 

W = ( ex P ( l^gSft ) for | sin 6| > 0.75; (1) 
I- otherwise 

which vanishes smoothly at sin \ b\ = 0.75 or |6| = 48.6°. Thus only the high latitude 25% of sky is 
used. 

A 2-D fast Fourier transform is then applied to this windowed map. Since the 512 pixel 
size of the image covered a hemisphere or 180°, the frequency spacing in the output nominally 
corresponds to A£ = 2. However, the use of an equal area projection instead of an azimuthal 
equidistant projection expands the central portion of the map by a factor tt/V8, so the actual 
output step corresponds to A£ = 2v / 8/tt = 1.8. The power spectra of the Northern and Southern 
polar caps were added and then summed into rings of constant |fc|. Finally the power spectrum 
of the difference map was subtracted from the power spectrum of the sum map to give the power 
spectrum of the sky. Figure || shows the resulting power spectra for DIRBE bands 7-10: 60, 100, 
140 and 240 /im. Points with £ < 19 come from the Wright et al. (1994) Hauser- Peebles method 
applied to the sky with \b\ > 48.6°, while the FFT method gives points with £ > 9. The agreement 
in the overlap region is good. 



5. Comparison with IRAS 

The power spectra of several cirrus clouds observed by IRAS are given by Gautier et al. 
(1992) in the form 

p(p) = p(o)(p/ Po r (2) 

with p Q = 0.01 cycles/arcmin and P(0) ranging from 4 x 10 7 to 4 x 10 10 Jy 2 /sr. In order to 
convert this to CV's I first note that p Q corresponds to £ = 216 since £ is the number of cycles in 
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Fig. 4. — Power spectra of the DIRBE 60, 100, 140 and 240 fan maps in |6| > 49°. The noise 
spectra for 140 and 240 pm are flat, as expected for a detector noise dominated map. The noise 
spectra for 60 and 100 pm have a pronouced shoulder at £ = 26, which is due to residual zodiacal 
light modeling error in alternate weeks. 



2n radians = 21600'. The variance of the sky in a band of frequencies of width dp = p centered 
on p is 27rpoP(0) = 7425 sr _1 P(0). The same variance computed from Ci is (Air)~ l £{2£ + l)Ci at 
£ = 216. These variances are identical so 

C i = ^\p/lf[l + (2£)- l ]P(p). (3) 

The factor (p/£) 2 = (47r 2 sr) _1 so the net result is that Q in (MJy/sr) 2 is numerically equal to 
P{p) in MJy 2 /sr except for the negligible [1 + (2£) -1 ] correction factor. 

I estimated the values of P (0) at p Q = 0.01 cycles/arcmin for the power spectra shown in Figure 
| using P(0) = median{(£/216) 3 C £ }. The values obtained are P(0) = 2.2 x 10 5 , 1.0 x 10 6 , 3.6 x 10 6 , 
and 2.1 x 10 6 Jy 2 /sr at 60,100, 140, and 240 pm. The value at 100 pm is considerably below the 
range seen by Gautier et al. (1992) in cirrus clouds, but it is higher than their extrapolation to 
dim cirrus. Since most of the sky with \b\ > 48.6° is quite dim, the IRAS and DIRBE power 
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Fig. 5. — Power spectrum amplitude C21& vs. wavelength (points) compared to the square of the 
galactic flux. 

spectra agree reasonably well, with the DIRBE measurement for |6| > 48.6° only 2.5 x the Gautier 
et al. (1992) estimate for |6| > 50°. 

If the cirrus clouds have the same geometry at all wavelengths, then the power spectrum 
should scale with frequency like C216 oc P(0) oc F% where F u is the flux of a typical cloud. Using 
the flux of the Milky Way for F v , I find that 

C216 = P(0) = W' 9 - u Fl (4) 

This scaling is shown in Figure ||. 

The observed power spectrum can be used to estimate the typical cirrus contribution to the 
noise using the tables in Gautier et al. (1992). With an aperture of 1.2X/D, and subtracting 
two reference apertures each two aperture diameters away, the expected cirrus noise at 100 fj,m is 
0.23 mJy for SIRTF and 0.54 mJy for ISO. These estimates are 40% smaller than the estimates 
in Gautier et al. (1992). My calculation requires a large extrapolation in angular scale from the 
0.7° DIRBE beam to the 0.5' SIRTF beam, but no extrapolation to the typical high latitude sky. 
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The cirrus noise at 240 /mi in a 3 meter submillimeter telescope such the planned FIRST mission 
using the same double subtraction scheme is 0.11 mJy. The Milky Way would give a flux 8 times 
greater than this cirrus noise at z = 0.7 for H Q = 50 km/sec/Mpc and = 1. 



6. Discussion 

The angular power spectra of the high galactic latitude infrared sky observed by DIRBE 
show the same power law behavior seen in the IRAS maps: Ce oc This is one power of I 
steeper than the power spectrum expected for the cosmic microwave background, so the problem 
of galactic contamination in moderate resolution experiments will be smaller than the galactic 
contamination in the COBE DMR maps. Thus proposed new CMB anisotropy experiments that 
observe with ~ 0.5° beams will not be ruined by contamination from galactic dust. 



The DIRBE data products are the result of dedicated work by hundreds of scientists and 
engineers who made the COBE project a success. I have had many useful discussions with W. 
Reach, T. Kelsall, M. Hauser and other DIRBE team members about zodiacal light modeling. 
However, the particular zodiacal light model used in this paper and especially its errors shown in 
Figure [l] are entirely my own. 



A. Appendix: Zodiacal Light Modeling 



The models used in this paper are examples of physical models that integrate a volume 
emissivity along the the line of sight specified by unit vector n: 



h{n,t)= p{r{s)) KiDiiTiR)) + ai^)Di(TQ)R- 



ds 



(Al) 



where R = \r\, Ki is a coefficient giving the efficiency in the i th DIRBE band of the interplanetary 
dust grains for thermal emission, while gives their scattering efficiency. &(p) is the phase 
function for scattering in terms of /i, the cosine of the scattering angle. Di(T) is the integral of the 
Planck function through the DIRBE filter response functions. The position vector r is given by 



r(s) = a?©(i) + ns 



(A2) 



where XQ)(t) is the position of the Earth in its orbit around the Sun. The density p(f) and 
temperature T(R) of the interplanetary dust cloud, and the scattering phase function ^(/i) 
are all represented by functions with adjustable parameters pi having values listed in Table |]. 
Scaled values of the grain efficiency factors Ki and Oi are listed in Table || with scalings given by 
Equations A17 and A18. The values of all the parameters are determined by assuming that all of 



the observed time variation in the DIRBE signal for a given pixel n is due to the motion of the 



Earth around the Sun described in Equation A2. 
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However, the symmetry axis of the zodiacal cloud is not the ecliptic pole but is both tilted 
and offset (Dermott et al. 1994). The tilt is parameterized by the pole unit vector given by 

zc = / P6 ' P7 ' 10) (A3) 

The inclination i with respect to the cloud is specified by 

sini = ^-^. (A4) 
R 

The density contours in the cloud are also not centered on the Sun. An offset radius is defined as 

R c = R + o c -f (A5) 

where 

o c = (p 8 /10,p9/10,0). (A6) 
The simplest models for the zodiacal cloud assume 

p{r) = Po f(i)R- pi (A7) 

where f(i) specifies the heliocentric inclination dependence of a fan shaped cloud and p\ is the 
exponent of a power law density dependence. The density scaling p Q is fixed at unity and the 
overall strength of the zodiacal signal is determined by the «j's and the <7j's. The temperature 
dependence is also given by a power law: 

T(r) = T R- p ' 2 (A8) 

The temperature was fixed at R = 1 au at 280 K because of a strong degeneracy between the k^s 
and T a . The inclination dependence f(i) is given by 

f(i) = exp(— P3 x Z PA ) + P15 + piQ x sin 2 i + pn x 4sin 2 i exp(— 4sin 2 i) 

+ pis x 16 sin 2 i exp(— 16 sin 2 i) + pig x 36 sin 2 i exp(— 36 sin 2 i)) (A9) 



where Z is "rounded" absolute value of sini: 

^ j | sin ^ | — for | sini| > S; 

The break point is given by S = exp(ps). 



1 sin 2 i/S' for|sini|<5. 



(A10) 



The ring of particles resonantly trapped by gravitational perturbations from the Earth 
(Dermott et al. 1994; Reach et al. 1995) is modeled using coordinates in a frame rotating to follow 
, the mean orbital longitude of the Earth: 

xd = xcosL0 + y sinL0 
y D = y cos L0 — x sin 

ZD = z (All) 
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Parameter 


Value 


Description 


7)1 

p 1 


1.2186 


radial Hpnsitv pxnonpnt, 


P2 


0.4451 


radial t,piT)r)pra,tiirp pxnonpnt 


Pi 


3.6122 


vprtiral "spalp hpifrit" 


PA 


D 9285 


vprtipal npvimtv pvooriPTit 




-1.4766 


ln(sin i) at break 


P% 


0.3705 


10 x cloud pole x component 


Pi 


-0.0736 


10 x cloud pole y component 


PS 


-0.0235 


10 x cloud offset x component 


P9 


-0.0081 


10 x cloud offset y component 


Pw 


0.7548 


10 x density contrast of Dermott ring 


Pll 


5.5301 


ln(Te), band temperature at R = 1 


P12 


0.7727 


scale factor for bands 


P13 


0.4284 


the "dimple" in Dermott's ring 


PU 


27.7741 


vertical scale for Dermott's ring 


P15 


-0.0251 


spherical term in vertical density 


PIG 


0.0249 


sin 2 i term in vertical density 


Pn 


-0.0456 


additional density at sini rs 0.5 


P18 


-0.1276 


additional density at sini ~ 0.25 


P19 


-0.0103 


additional density at sini ~ 0.17 


P20 


-0.3133 


phase function linear coefficient 


P21 


0.5749 


phase function quadratic coefficient 



Table 1: Diffuse Cloud Parameters 
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Then the following terms are defined: 



Ld = |atan2(yo, xp) + 0.25 1 

cos(87tLd/3) for < |; 



A 



D 



(cos(87rL D /3) - l)/2 for § < L D < §; 







exp 



otherwise. 



56.5 Jxf, +y 2 D - 1.133 + 0.133 Pl3 exp(-4L£) 



2 ^ 2 Pl4^D 



i? 2 



(A12) 
(A13) 

(A14) 



The final form for the density is 



^) = ^R P1 (1 + ^ 



The phase function for scattering is given by 

= exp(-p 2 o/^ + P2iM 2 ) 



(A15) 



(A16) 



The values of Ki and a% are also adjustable parameters. Only bands 1-3 (J, -fT and L) have 
non-zero <7j's, with 



Oi = 1.8p 2 i+i x 10 
and only bands 3-10 have non-zero k^s with 

«i = 10 _7 p22+i- 
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(A17) 
(A18) 



However, this model does not include the IRAS bands (Low et al. 1984) The bands are 
described by more parametrized functions with parameter values given in Table [3|. So the following 
term is added to the intensity: 



with 



Ib=P12 I PB{r{s)) [^(TBiR)) + oM^)Di(Tq)R- 2 
T B (R) = - 6XP(P11) 



ds 



R 



The offset and tilt of the bands are given by 

(95)96,10) 



Z B 



2 i „2 



100 + q£ + 



(A19) 
(A20) 

(A21) 



where q n will denote adjustable parameters in the models. The inclination ig with respect to the 
bands is specified by 

zb • r 



8111 1 B 



R 



(A22) 
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Table 2. Scattering and Emission Efficiencies 



Parameter Value Description 



P22 


0.6864 


0"i scaling 


P23 


0.8362 


0"2 scaling 


£>24 


0.6346 


o"3 scaling 


P25 


2.2829 


K3 scaling 


P26 


1.4175 


«4 scaling 


P27 


1.2422 


K5 scaling 


P28 


1.2996 


kq scaling 


P29 


0.8788 


k-j scaling 


P30 


0.7646 


Ks scaling 


P31 


0.8243 


Kg scaling 


P32 


0.6254 


Kio scaling 



Table 3. IRAS Band Parameters 



Parameter 


Value 


Description 


qi 


1.3849 


10 x (smi) max for Band 1 


Q2 


0.2743 


Band 1 normalization 


Q3 


0.2807 


10 x (smi) max for Band 2 


Q4 


0.4407 


Band 2 normalization 


Q5 


0.1735 


10 x band pole x component 


Q6 


-0.2088 


10 x band pole y component 


Q7 


-1.5723 


10 x band offset x component 


Q8 


-0.2225 


10 x band offset y component 


Ri 


3.14 


Outer radius for Band 1 


R2 


3.02 


Outer radius for Band 2 
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The density contours in the bands are also not centered on the Sun. The offset radius is defined as 

R B = R + o B -f (A23) 

where 



o B = (97/10,98/10,0). 



The density in the bands is given by 

2 



R V- [fcoshfi^i 
-n-B 7= i L L 



sin i£ 

T 



for | sini^l < q2j-i and -Rb < Rj\ 
otherwise. 



(A24) 



(A25) 



The functional form cosh(1.72x) was chosen instead of (1 — x 2 ) -1 / 2 to avoid numerical 
problems when integrating through the cusps. Note that the coefficient is chosen so that 
j^l cosh(1.72x)<ix = it = f^(l — x 2 )~ l l 2 dx. The outer radii are taken from Jones & 
Rowan-Robinson (1993). 

The band parameters {qj} and the diffuse cloud parameters {p{\ are optimized in alternating 
runs. When optimizing the band parameters a set of normal points emphasizing the ecliptic is 
used. When optimizing the diffuse cloud parameters a separate set of normal points emphasizing 
the ecliptic poles is used. 

A second model that used a slightly different treatment of the Dermott ring and clump was 
also tried. An additional parameter Pi describes the radial extent of the ring: 



D = exp 



-56.5Pi [\Jx 2 D + y 2 D - 1.133 + 0.133pi 3 exp(-4L^; 



2 _ PUZ 2 D 

R 2 



(A26) 



A second additional parameter P2 gives a separate adjustment for the density of the "clump" 
trailing the Earth. The density of the cloud is given by 



P(r) = ^-jmc Pl (l + 



Pl0 D + P 2 DA) \ 
10 ) 



(A27) 



Fitting for these parameters gives P2 > p±o which implies a negative density in the ring+clump 
system, which seemed unphysical, so the simpler model described above was developed. 
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